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PREFACE 


This  Memorandum  has  been  prepared  as  a  part  of  a  continuing  study 
of  hypersonic  lifting-vehicle  technology.  It  is  an  extension  of  pre¬ 
vious  investigations  of  the  use  of  aerodynamic  forces  to  provide 
lateral-range  capability  for  glide-reentry  vehicles.  The  closed- 
form  solutions  presented  here  should  be  of  interest  to  designers  and 
planners  concerned  with  the  preliminary  design  and  capabilities  of 
hypersonic  glide  vehicles. 
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SUMMARY 

Approximate  closed-form  solutions  for  various  flight  conditions 
have  been  obtained  to  determine  both  the  longitudinal  and  lateral  range 
of  hypersonic  glide  vehicles  with  constant  bank  angle.  Results  for 
equilibrium-glide  vehicles  with  a  constant  lift-to-drag  (L/D)  ratio 
and  small  and  slowly  changing  flight-path  angle  are  presented  in  graph¬ 
ical  form.  Other  approximate  closed- form  solutions  are  also  obtained 
for  glide  reentry  at  very  small  flight-path  angle,  near-constant- speed 
glide  at  high  altitude,  constant-deceleration  glide  at  constant  alti¬ 
tude,  and  constant-deceleration  glide  at  fixed  flight- path  angle. 

The  assumpMon  of  a  very  small  flight-path  angle  (yw  0)  results 
in  a  smaller  range  prediction  than  does  the  assumption  of  a  small  and 
slowly  changing  flight-path  angle  (y  «  0.01  rad).  This  is  especially 
true  for  prediction  of  the  lateral  range  of  vehicles  with  high  lift- 
to-drag  ratio.  For  a  vehicle  with  an  L/D  of  3,  entering  at  0.98  or¬ 
bital  velocity  and  decelerating  to  0.2  orbital  velocity,  the  assump¬ 
tion  of  a  very  small  flight-path  angle  results  in  predictions  that 
underestimate  lateral  and  longitudinal  range  by  22.6  and  22.5  percent, 
respectively,  in  comparison  with  the  predictions  based  on  the  assump¬ 
tion  of  a  small  and  slowly  changing  flight-path  angle.  In  other  words, 
completely  neglecting  the  flight-path  angle  in  the  equations  of  motion 
leads  to  a  conservative  range  prediction. 
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SYMBOLS 


A  *  area 


B1  =  Dl(cos  §/Q 


C  =  constant  o£  integration 
Cp  =  drag  coefficient 

C_  =  lift  coefficient 
L 

2 

D  =  drag  force  =  C^pu  A/2 

dd  -  Vo'^V) 

CL  - 

g  =>  acceleration  due  to  gravity 
h  =*  altitude 
he  *  entry  altitude 

h.  **  initial  altitude 

i. 

L  *  lift  force  *  C^pu^A/2 
Rg  =  earth's  radius 


r  •  dimensionless  altitude,  Eqs.  (7) 

%  ’  V(*E  +  V  “  W 
*1  *  V<RE  +  V  "»  hi/SE 

S  ■  distance  along  the  flight  path 
t  *  time 

u  ■  flight  speed  at  any  given  time 
uft  »  entry  speed 

«  initial  speed 


‘■JlSifrfHfffl 
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cu  » 


orbital  speed  at  reference  altitude  or  entry 

velocity  ratio  «  u/uo>  Eqs.  (7) 

reentry  velocity  ratio  *  u  /u 

e  o 

initial  velocity  ratio  *  u  /u 

i  o 

vehicle  weight 
downrange  distance 
aider ange  distance 
tan  #/[yO  -  2)] 

2/[y(1*/D)  cos 
flight* path  angle 

aean  value  of  small  ami  slowly  changing  flight- path  angle 
flight-path  angle  at  entry 

initial  flight-path  angle 

inverse  scale  height 

C(RE  +  he}  %  V 
lateral-range  angle 
atmospheric  density 
reference  atmospheric  density 

dimensionless  density  ratio,  Eqs.  (7) 

hank  angle 

turn  angle 

turn  angle  at  entry 

initial  turn  angle 
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i.  nmuffiocnoM 


The  problem  of  obtaining  closed-fora  solutions  for  glide-reentry 
vehicles  has  been  an  area  of  interest  for  many  years*  Since  the  general 
exact  analytical  solutions  of  the  fundamental  equations  of  notion  for 
glide  reentry  are  difficult  to  obtain,  various  approximations  have  been 
made  by  different  authors  for  different  flight  conditions*  For  instance, 
the  closed  solutions  for  very  small  angles  of  inclination  have  been  de- 

(1  2) 

rived  by  Gazley.  Eggers  et  al.  also  obtained  approximate  closed- 

(3) 

form  solutions  for  small  angles  of  inclination.  Chapman  has  devel¬ 
oped  analytical  solutions  for  small  angles  of  inclination  and  small 
lift- to- drag  (L/D)  ratios.'  J  The  solutions  for  moderate  angles  of 
inclination  and  medium  L/D  ratios  were  solved  by  Lees  et  al. Loh 
has  derived  closed- fora  solutions  for  negative  lift  at  small  and  large 
angles  of  inclination.^  For  entry  at  small  L/D  ratios  and  large  an¬ 
gles  of  inclination,  the  approximate  results  of  Arthur  are  available. ^ 
Loh  later  obtained  a  second-order  approximate  solution  for  a  noncscil- 
latory  type  of  entry  trajectory  and  derived  some  exact  analytical  solu- 
tions  for  certain  trajectory  variations.  *  J  Recently,  Cohen  has  de¬ 
veloped  the  closed- form  solutions  for  a  constant-deceleration  flight 
path.  Loh  extended  his  second-order  solutions  to  an  oscillatory 

type  of  entry  trajectory.  However,  most  of  these  results  are  for 

two-dimensional  trajectories  where  the  bank  angle  is  zero. 

The  closed-form  solutions  for  c  three-dimensional  trajectory  at  a 

constant  and  very  small  angle  of  inclination  with  constant  bank-angle 

(12) 

control  have  been  obtained  by  Inland.  This  Memorandum  extends 

Hyland* s  results  to  obtain  exact  closed- form  solutions  for  equilibrium- 
glide  path  at  small  and  slowly  changing  angles  of  inclination,  as  well 
as  for  some  special  cases  where  the  angles  of  inclination  are  larger. 

A  constant  bank  angle  is  assumed  for  both  flight  conditions.  The  re¬ 
sults  for  an  equilibrium-glide  path  at  small  and  slowly  changing  angles 
of  inclination  are  also  shown  in  graphical  form  to  provide  information 
for  preliminary-design  purposes. 
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The  general  equations  of  motion  with  constant  bank  angle  in  a 
nenrotating  three-dimensional  inertial  coordinate  system,  as  shown  in 
Fig.  1,  are 


~  u  -  sin  y  -  §  (1) 


"  g  ^  "  W  COS  $  ^  "  ^2^cos  Y  (2) 


a  a)  cos  Y  -  |  sin  *  -  sLsgg.2.  r  (3) 

®  u“ 

o 

where  u  is  the  speed  at  any  given  time,  t  is  the  time,  Y  is  the  flight- 
path  angle  in  the  vertica.  direction,  D  is  the  drag  force,  W  is  the 
vehicle  weight,  L  is  the  lift  force,  $  is  the  bank  angle,  uq  is  the 
orbital  speed  at  an  entry  altitude  of  400,000  ft,  cu  is  the  vehicle 
turn  angle  (heading  angle),  and  X  is  the  lateral-range  angle.  Equa¬ 
tion  (1)  pertains  to  motion  along  the  flight  path,  Eq.  C2)  to  motion 
up  and  down  the  local  vertical,  and  Eq.  (3)  to  a  direction  out  of  the 
flight  plane.  It  is  assumed  that  the  altitude  region  of  interest  is 
below  400,000  ft  and  that  gravity  does  not  change  very  much  between 
sea  level  and  400,000  ft.  It  is  further  assumed  that 

2  2 

L  .  *  s  u  cos  Y  cos  u)  tan  X 

D  2 

u 

o 


and 


sin  Y  «  Y  cos  Y  «  1 


Fig.  1— ' Three-dimensional  flight-trajectory  coordinates 
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The  first  assumption  indicates  that  the  centrifugal-force  component 
in  the  lateral  flight  plane  is  small  in  comparison  with  the  lift  compo¬ 
nent  in  the  lateral  flight  plane.  This  assumption,  which  is  used  through- 

(12) 

out  this  study,  is  also  used  by  Nyland.  However,  for  constant-alti¬ 

tude  glide  at  superorbital  speeds,  this  assumption  may  cause  an  error  of 

(13) 

30  percent  for  L/D  =  2  in  heading-angle  change,  as  pointed  out  by  Wang. 

The  second  assumption  has  been  relaxed  for  one  special  case,  given  later. 
With  these  assumptions,  Eqs.  (1)  to  (3)  can  be  written  as 


1 

S 


€ 

U 


Y 


(4) 


u  *  L 
—  v  ■  r:  cos  * 
g  w 


(5) 


u  •  L  .  * 

—  mar:  sin  * 
g  W 


Let 


dS 
*  dt 


dh  -  -YdS 


V 


h 

r  s*  ■  ■  ■  ~ 

*E  +  he 


(6) 


(7) 


where  S  is  the  distance  along  the  flight  path,  h  is  the  altitude,  he 

is  the  entry  altitude,  and  R_  is  the  earth's  radius.  These  variables 

“  -  (10) 

were  introduced  by  Cohen,  but  p  has  been  modified.  Equations  (4) 

to  (6)  can  then  be  rewritten  as 


dV2 
d  t 


-2 


(8) 
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GLIDE  REENTRY  TRAJECTORIES  AT  SHALL  AND  SLOWLY  CHANGING 
FLIGHT-PATH  ANGLE 

One  kind  of  atmospheric  entry  by  lifting  vehicles  at  small  angles 
of  inclination  is  an  equilibrium-glide  flight  path,  where  the  gravi¬ 
tational  force  is  balanced  by  the  lift  force  and  by  the  centrifugal 
force  due  to  the  curvature  of  the  flight  path.  With  the  assumption 
of  a  slowly  changing  flight-path  angle,  or  y  **  0  and  jyj  <  Y,  which  is 
also  assumed  by  Gazley,^  Eggers  et  al.,^  and  Nyland,  Eqs.  (4) 
to  (6)  reduce  to 


1 

g 


u 


Y 


(14) 


0  =*  —  cos  5-1 - -r  (15) 

w  i 

u 
o 


—  da  sx  ~  sin  5  (16) 

g  W 

Gazley^  and  Eggers  et  al.  ^  solved  Eqs  (14)  and  (J5)  for  the 

(12) 

two-dimensional  case,  where  §  =  0,  by  assuming  Y  aO;  while  Nyland' 
obtained  solutions  for  the  three-dimensional  case,  solving  Eqs.  (14) 
and  (15)  by  assuming  Y  «  0  and  5  =  constant.  Although  it  is  asst  ed 
here  that  5  =  constant,  Y  is  considered  in  Eq.  (14)  as  small  but  not 
negligible.  In  terms  of  the  new  variables  V  and  r,  as  introduced  in 
(7),  Eqs.  (14)  to  (16)  can  be  rewritten  as 


Y  + 


1  iid 

2  dr 


V2 


cos 


(14a) 


(15a) 

(16a) 
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Dividing  Eq.  (14a)  by  Eq.  (15a)  yields 


dV2 

dr 


+  3v2 


3  - 


2 


(17) 


where 


3 


2 _ 

cos  $ 


and  is  a  function  of  Y,  L/D,  and  $.  The  solution  of  Eq.  (17)  is 


-J » 

(3  -  2)  dr  +  Ce 


dr 


(18) 


where  C  is  a  constant  of  integration.  The  Initial  condition  i$  u  »  ug 

at  h  =  h  ,  or  V  =  V  at  r  =  r  .  Then  Eq.  (18)  can  be  rewritten  as 
C  c  c 


_  f  3  dr  -F  3  dr  f  3  dr 
r  »  f eJ-  +eJ  eJ  (3  -  2)  dr  (19) 

C  T 


For  constant  L/D  and  §,  and  small  and  slowly  changing  Y»  3  is  approxi¬ 
mately  a  constant,  or  a  mean  value  is  assumed  for  Y»  ari  Eq.  (19)  can 
be  reduced  to 

Equation  (20)  gives  the  relation  between  altitude  and  velocity  for  con¬ 
stant  values  of  3<  Since  3  is  a  function  of  flight-path  angle,  L/D  ra¬ 
tio,  and  bank  angle,  any  changes  in  one  of  these  quantities  will  affect 
the  vehicle  flight-speed  and  altitude  relationship.  Dividing  Eq.  (16) 
by  Eq.  (15)  and  taking  3  and  y  at  mean  values,  respectively,  one  can 
obtain  the  turn  angle: 
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U) 


tan  f 

9(8  "  2) 


2 

F 


in 


8V2  -3  +  2 

SV2  -  3  +  2 
e 


+  in-§ 
V 


(21) 


or 


(Ji 


tan  $ 
Y(3  -  2) 


2(re  -  r) 


-  in 


-I 


lii-j-j) 


3(Vr) 


(22) 


In  the  flight  regime  of  interest,  Y  is  small  and  slowly  changing. 
Therefore  it  is  assumed  that  a  mean  value,  y,  can  be  obtained.  The 
initial  turn  angle  is  taken  as  zero  in  Eqs.  (21)  and  (22).  Dividing 
Eq.  (14)  by  Eq.  (15).  one  obtains 


1  .  _ 

gU~  2  L 


u  «  Y  - 


U  ~  C08  $ 

o  D 


L  £ 
D  c°‘  * 


(23) 


Solving  Eq.  (23)  for  t  as  a  function  of  V  by  assuming  constant  L/D 
and  •  and  mean  value  Y. 


"o  d  co* 1  rir  [v^-v  s  -  2][ve# +ve  -  2] 

2g  Vs  -  2  in  LW6  +V8  -  2][veVe  -  V6  -  2]  <24) 


Hoting  that 


.  1  du2 

u  *  2  dS 

where  S  is  the  distance  along  the  flight  path,  Eq.  (23)  can  be  re¬ 
written  as 


1  du2 
2g  dS 


2  L 

u  —  cos 
o  D 


$ 


2 

u 


Y  - 


L  . 
D  co*  * 


(23a) 


-9- 


Solving  Eq*  (23a)  for  S  as  a  function  of  V  and  r,  respectively,  one 
obtains 


Ro  .  V2^  +  2-8 

-  —  in  ~7r - tL 

BY  7*8  +  2  - 


and 


S 


-  r) 


(26) 


wbere  R 

o 


2. 

“o/g- 


The  downrange  distance  can  be  obtained  from 


dx  =  dS  cos  u: 


(27) 


where  x  is  the  distance  traveled  in  the  original  direction  of  motion 
and  can  be  taken  as  approximately  the  distance  on  earth  because  of  the 
small  flight-path  angle  and  small  ratio  of  altitude  to  earth  radius. 
For  constant  §  and  L/D  and  mean  value  Y,  Eq.  (27)  can  be  integrated  as 


-R  r 

X  =  r  COS  •( 

Y  J- 

r 

e 


2(re  -  r)  -  xn 


.  <  4-i-  4  -’6n 


>dr 


(28) 


In  terms  of  7 


cos  y  or 


R  V 
o  - 

x  =  ~  r 

V  I 


2  ,  3V2  -  3  +  2  ^ 

a  xn  a  xn  . 

p  3\r  -  3  -i-  2  v 


v  3  +  2  -  3 


dV 


(29) 
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Simi larly  the  siderange  distance  can  be  obtained  for  constant  $  and 
L/D  and  mean  value  V  from 


dy  ■  dS  sin  uj 


(30) 


where  y  is  the  distance  traveled  perpendicular  to  the  original  direc¬ 
tion  of  motion,  or 


y 


sxn 


a 


r) 


£n 


In  terms  of  V 


fJ 


.V2 


sin  y 


or 


#iafV*-B  +  2  +  x, 


8  ev2  - 


B  +  2 


V2 

e 

V2 


V2 

e 


+  2  -  £ 


dV2 


(32) 


where  or  =  tan  §/(&  -  2)v 


GLIDE  REENTRY  TRAJECTORIES  AT  VERY  SMALL  FLIGHT-PATH  ANGLE 

For  equilibrium  glide  at  very  small  flight-path  angle,  such  as 
for  high-L/D  flight,  Y  «  Y  0.  Equations  (4)  to  (6)  can  be  reduced 
to 


u  _  D 
g  ~  » 


w  cos  * 


U 

O 


(33) 


(34) 


_n  _ 


u  •  L  ,  . 

—  to  *  r:  sin  f 
g  W 


(35) 


From  the  definition  of  3 


■§•  =  Y  “  cos  v  =  0 


(36) 


Substituting  Eq.  (36)  into  Eqs.  (21),  (24),  (25),  (29),  and  (32), 
respectively,  one  obtains 


L  .  *  .  i 

*  =  d  Sin  *  in  V" 


(37) 


u  —  cos  §  (V  -  1)(V.  +  1) 
t  -  -2-2= -  in - 1 - 


2g 


(V  +  1)(V.  -  1) 


(38) 


S  *  |  R  |  cos  §  in  1  ~ 

2  oD  1  +  V? 


(39) 


cos 


1  _  L 

X  »  -T  R  —  cos  * 

2  o  D 


■r 


V?' 


1L  .  ~ 

2Dn”U7 


dV2 


V2- 


(40) 


iB  l  .  r 

o  R  n  cos  * 

£  o  u 


r  2“ 

v~ 

dV2 

sin 

—  ^  sin  $  in  —x 

L2  D 

'v2 


V2- 


(41) 


where  V  * 


u  /u  ,  and  u.  =  initial  speed, 
i  o*  x 
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Equations  (33)  and  (34)  and  their  results,  Eqs.  (37)  to  (41),  are  the 

(12) 

same  as  obtained  by  Kyland. 

NEAR-COKSTAKT-SPEED  GLIDE  AT  HIGH  ALTITUDE 

For  some  applications  of  atmospheric  flight,  it  may  be  desirable 
to  have  the  lifting  vehicle  enter  at  near-orbital  speed  and  fly  at 
maximum  L/D  ratio  or  maximum  C  and  at  constant  bank  angle  until  it 

la 

reaches  denser  atmosphere  at  an  altitude  of  about  250,000  ft.  In  this 
phase  of  flight  at  high  altitude,  atmospheric  drag  is  assumed  to  be 
nearly  equal  and  opposite  to  the  component  of  gravity  along  the  flight 
path.  Thus,  there  would  be  little  change  in  vehicle  speed  and,  in  the 
present  notation,  this  implies  that  l/V2-!  «  0.  Let  the  atmospheric 
density  be  expressed  by  the  well-known  exponential  approximation  or 


P  - 


(42) 


where  pQ  is  the  reference  atmospheric  density  and  £  is  the  inverse 
scale  height.  Substituting  Eq.  (42)  into  Eq.  (7), 


P 


poUo  -£h 
W/A  C 


(43) 


where  £  =  £(Rg  +  hg)  «  £8^ 

Then  for  flight  at  small  angles  of  inclination,  Eqs.  (4)  to  (6) 
can  be  reduced  to 


u 

g 


Y 


u  -  L  _ 
-  -  Y  -  -  co.  « 


u 

I 


03 


wsin* 


(44) 

(45) 


(46) 
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Or  in  terms  of  V  and  r,  Eqs.  (44)  to  (46)  become 


2 

dV2  1.  PoUo  -Qr  yZ  2 

df  Y  W/CdA  e 


(44a) 


2 

PoUo  -Cr 

w/cla 


cos  x 


(45a) 


du? 

df 


i  p„uo 

2Y  W/CjA  e 


sin  § 


(46a) 


Let 


P  u2 
o  o 

W/CjA 


2 

Pouo 

w/cla 


Then  Eqs-  (44a)  to  (46a)  can  be  rewritten  as 


dV2 

df 


dti‘ 

df 


-  D  e'^1  V2  -  2 
Y  D 


_  -Cr  , 
D,e  cos  $ 


2Y  V  sin  * 


(44b) 


(45b) 


(46b) 


For  constant  W/CLA,  Eq.  (45a)  can  be  readily  integrated  to  yield 

-Lv 


14- 


Y2  +  d.  £22-1  L 

*  i.  l  L- 


-  -Cr 


where 


Vg  -  flight-path  angle  at  entry 
~  dimensionless  altitude  at  entry 

2  "Cr 
A1  "  Ye  +  V 

*  Dl(cos  i/Q 


Hence 


*  -  tv  v*t 


TAi  -  ^ 


1  *  A. 


1  ~  f  aT  e'Cr  +  •• 


'1  -Cr 


e  =  «  1 


Mr  -  ♦  v**) 4  -  £ 


Substituting  Eq.  (48)  into  Eqs.  (44a)  and  (46a),  respectively, 

dV2  1  ? 


2 
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— r  +  ~7-='  D  e  °r  sin  $  ~  0 
dr  27A  L 


Integrating  Eqs.  (49)  and  (50)  for  constant  W/C^A, 


V2 


V2e 

e 


1  (  “fr, 
c 


-r-t 


-  e 


D  1  e“£r  -  D  1  e“£r 

VTX  c  re  V3q  ? 

+  2e  e  ^ 

Jr 


u) 


£ 


L  sin  5  j^e  -  e  ”  eJ  +  a). 


(50) 


(51) 


(52) 


where  and  co^  are  initial  values  of  reentry  velocity  ratio  and  re¬ 
entry  turn  angle,  respectively. 

CONSTANT-DECELERATION  GLIDE  AT  CONSTANT  ALTITUDE 

After  rapid  and  close  to  constant  velocity  descent  from  orbit  to 
an  altitude  of  about  250,000  ft,  many  entries  of  interest  may  require 
the  vehicle  to  fly  at  constant  altitude  to  perform  plane-change  ma¬ 
neuvers.  In  order  to  maintain  vehicle  glide  at  constant  altitude,  the 
gravity  has  to  be  balanced  by  lift  force.  As  the  vehicle  decelerates, 
the  lift  coefficient  has  to  increase  to  overcome  the  reduction  in  dy¬ 
namic  pressure  until  the  lift  force  is  too  small  to  sustain  constant- 
altitude  flight.  In  this  phase  of  flight  h  is  essentially  constant 
and  y  *  Y  *  0*  The  governing  equations  of  motion  are  identical  to 
that  for  equilibrium  glide  at  very  small  flight-path  angle  except  that 
L  and  D  are  variables,  namely. 


u 

g 


(53) 


tallfr'l'li 
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L 

W 


/  i»  i  V 


u  ■  I*  .  - 

—  a)  =  —  sin  9 
g  W 


(55) 


However,  in  order  to  prevent  skip,  the  vehicle  is  required  to  fly  at 
variable  L/D  at  almost  constant  altitude.  In  addition,  the  vehicle 
may  fly  at  certain  fixed  bank  angles  to  achieve  desired  plane  changes. 
For  manned  maneuverable  vehicles,  it  may  be  desirable  to  perform  this 
aerodynamic  maneuvering  in  a  constant-deceleration  mode  of  flight. 
Then  u  and  D  become  constants.  Therefore,  it  is  assumed  that  the  L/D 
ratio  follows  the  relationship  suggested  by  Cohen,  namely. 


Substituting  Eq.  (56)  into  Eq.  (54), 


(56) 


W 

D  *  A2  cos  $ 


(57) 


or 


2  - 

V  A^P  cos 


(58) 


where  A2  *  -1/u/g  cos  $  is  constant  and  deceleration  (-u)  can  be  spe¬ 
cified  in  g*s.  Equation  (56)  can  also  be  rewritten  as 


Then  Eq,.  (53)  can  be  rewritten  as 


(59) 
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V  -  —  =  constant 
u 

o 


Therefore  flight  is  at  a  constant  rate  of  deceleration,  and 


V  = 


v.  +  £■ 

i  u 


t 


Solving  Eq.  (55),  one  obtains 


u> 


.  g  tan  § 

<i».  +  a - 

x  u 


t2  In  V  +  -  (V.  +  V) 
u  2  x 


(60) 


(61) 


or 


a»  = 

l 


fr  tan? 


(62) 


For  a  bank  angle  $  -*  90  deg,  Eqs.  (61)  and  (62)  approach  infinity  and 
the  results  are  not  applicable,  since  equations  of  motion  are  for 
$  <  90  deg.  The  vehicle  range  can  also  be  obtained  from  Eq.  (53); 


2 

s  ■  is  ^ 


(63) 


CONSTANT -DECELERATION  GLIDE  AT  FIXED  FLIGHT-PATH  ANGLE 

For  manned  reentry  flight,  it  may  be  desirable  to  maintain  con¬ 
stant  deceleration  during  the  high-speed  portion  of  flight  after  con¬ 
stant-altitude  glide.  If  the  flight-path  angle  is  fixed,  then  Eqs.  (1) 
to  (3)  can  be  rewritten  as 

—  *  sin  v  -  77  ■  constant  (64) 

e  w 
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0  = 


L 


w  cos 


^cos 


V 


(65) 


u  «  L  ,  . 

—  a)  cos  y  =  —  sin  * 

g  W 


(66) 


Equations  (64)  to  (66)  can  be  rewritten  in  terms  of  new  variables 
V  and  r  as 


dvf  V  2  sin  y 

dr  “  Y  V  "  Y 


(64a) 


0 


i  C.  PV2 


2  L 


CQ8  $ 

cos  Y 


-  1 


-v2 


(65a) 


dm  sin  $ 

dr  “  “  2Y  cos  Y 


(66a) 


If  the  same  L/D  ratio  expression  suggested  by  Cohen^0^  is  assumed, 
namely. 


(56) 


substituting  Eq.  (56)  into  Eq.  (65)  gives 


^  W  cos  Y 

D  »  - - 1— 

A.  cos  t 


(67) 


or 


,  ^  2  cos  Y 

®  V^A^P  cos  $ 


(68) 
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A2  "  /  if 


-cos  Y _ 

sin  cos  $ 


constant 


Equation  (66)  can  also  be  rewritten  as 


2A0  cos  Y  / .  .2 

c>  -  — r  H1- 

L  p  COS  9  \ 


-2  cos 


;(f- 


p[—  -  sin  Y  )  cos  $ 


2  $  V  V2 


Substituting  Eq.  (67)  into  Eq.  (64)  and  integrating  with  respect  to  t. 


V  =  V.  +  —  t 
i  u 

o 


Solving  Eq-  (64a)  for  V  as  a  function  of  r  gives 


2 

p  u  C_A  --  | 

»  oo  D  -£r  -  - 

2  sin  V  J  ^  W  T  1  'I 

Y  J-  6 

r 


P  u2  C  A 
,  o  o  D  -Cr 

Y  W  e 


jr  PoUo  SA  -Cr  - 

f  —IT  6 


+  V7  e 
i 


Dividing  Eq.  (64)  by  Eq.  (66),  one  obtains 


-  in.  +  S1°  i  $  2n  i  +  V?  -  V2 

1  2  “L  V?  1  J 


Equation  (64)  can  be  rewritten  as 
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III.  DISCUSSION  AND  CONCLUSIONS 


The  general  closed-form  solutions  presented  in  Eqs.  (11)  through 

(13)  can  be  solved  numerically  if  the  quantities  W/C  A,  H/C  A,  p,  and 

u  JLi 

as  functions  of  altitude,  are  known.  For  special  cases,  such  as 
equilibrium-glide  trajectories  with  small  and  slowly  changing  flight- 
path  angles,  the  closed -form  solutions  for  constant  L/D  ratio  and 
bank  angle  are  given  in  Eqs.  (20)  through  (22),  (24)  through  (26),  (28), 
(29),  (31),  and  (32)  in  terms  of  the  parameter  8  =  2/[y(L/D)  cos  §]. 

One  can  directly  obtain  the  value  of  8  for  various  combinations  of  Y, 
L/D,  and  cos  $  from  Fig-  2.  When  Y  and  §  are  fixed,  8  decreases  with 
increasing  L/D.  (The  range  of  8  which  is  of  interest  varies  with  L/D 
as  shown  in  Fig.  9.) 

All  the  results  shown  in  Figs.  3  to  9  are  for  equilibrium  glide, 
where  the  Y  term  is  retained  in  the  equations  of  motion  and  is  taken 
at  mean  value.  For  conditions  where  L/D,  Y,  and  4  are  not  constants, 
one  can  divide  the  flight  path  into  intervals  and  can  apply  all  the 
results  by  using  mean  values  of  L/D,  Y,  and  4  in  each  interval.  Total 
flight  time  and  total  range  can  be  obtained  by  summing  up  all  the  in¬ 
tervals. 

The  flight-speed  to  entry-speed  ratios,  as  functions  of  altitude 
ratios  for  different  values  of  8  as  given  by  Eq.  (20),  are  presented 
in  Fig.  3.  For  a  given  8  value,  the  speed  ratio  is  fixed  for  a  given 
altitude.  The  results  also  indicate  that  the  flight  speed  decreases 
with  decreasing  altitude.  However,  the  rate  of  decrease  is  slower  for 
the  higher  values  of  8.  For  given  valuer  of  L/D,  Vg,  and  small  and 
slowly  changing  Y,  the  trajectory  becomes  steeper  with  increasing  back 
angle.  A  high- L/D  vehicle  will,  in  general,  fly  a  steeper  equilibrium- 
glide  path. 

Figure  4  shows  the  influence  of  entry  speed,  bank  angle,  and 
flight-path  angle  on  turn  angle  a*  expressed  by  Eq.  (21).  The  re¬ 
sult  indicates  that  the  turning  rate  increases  as  the  vehicle  slows 
dcwn.  Large  turn  angles  can  be  accomplished  with  high  L/D  and  steep 
bank  angle.  The  effects  of  entry  speed  on  turn  angle  are  insignifi¬ 
cant  at  low  values  of  flight  speed  and  small  values  of  Y(L/D)  sin  4. 


Fig. 2— Relation  between  flight-path  angle,  L Id  ratio, 
bank  angle,  and  parameter  B 
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Fig. 3— Reentry  glide  path  for  constant  £ 
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Ve  =  0.98 
7  =  0.01  rod 
4>  =  30  deg 


Fig.5a — Turn  angle  and  altitude  relation 
(Vo  =  0.98,  /=  0.01  rad,  <£=30  deg) 
(Eq. (22)) 


Fig. 5b— Turn  angle  and  altitude  relation 
(Ve  a  0.98,  /  =  0.01  rad,  <$=45deq) 
(Eq.  (22)) 
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Fig. 5c— Turn  angle  and  altitude  relation 
(Ve  =0.98,  y  =  0.01  rad, $-60 deg  ) 
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Fig. 5d— Turn  angle  and  altitude  relation 
(Ve  =0.98,  r  =  0.03  rad,  $=45  deg) 


Fig. 5e — Turn  angle  and  altitude  relation 
(Ve  .  0.90,  /“  0.01  rad,  <$-  45  deg ) 
(Eq.(22)> 
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Fig. 7a— Velocity  and  flight-path -length  relation  (7=  0.01  rad) 
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Fig. 9a — Influence  of  bank  angle  on  ground  trace 
(L/D  =  1,  Ve  =  0.98) 

(Eqs.  (29),  (32)) 


Fig. 9b—  influence  of  bank  angle  on  ground  trace 
{L/D=2,  Ve  =  0.98) 

(Eqs  *29),  (32)) 


Fig.9c— Influence  of  bank  angle  on  ground  trace 
(L/D  =  3,  V8  =  G.98) 

(Eqs.  (29),  (32) ) 
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Fig.9d— influence  of  bank  angle  on  ground  trace 
( L/D  =  4,  Ve  *  0.98) 

(Eqs.{29),  (321) 


rig.9e—  influence  of  bank  angle  onground  trace 
( L/D  =3,  Ve  =  0.90) 

(Eqs.  (29),  (32)) 
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The  altitude  and  turn-angle  relation  as  obtained  from  Eq.  (22) 

i  *  T>i _  C-.  «.vw _ t-  C  -  A  _ _ 

to  OUUWU  XU  X'  J-  g  O  •  ja  WiUUU^U  JC*  AO  CA^C^VCU)  WUO  bux  U  au^XCi  canto 

with  increasing  bank  angle  and  the  turning  rate  is  highest  at  the  end 
of  a  turn  for  constant  bank  angle.  A  comparison  of  Figs.  5b  and  5d 
shows  that  for  given  flight  conditions,  the  turn  angle  is  larger  for  a 
small  flight-path  angle.  Figures  5b  and  5e  show  that  for  fixed  3,  bank 
angle,  flight-path  angle,  and  velocity  ratio,  the  lower- entry- velocity 
vehicle  will  achieve  a  given  turn  angle  at  higher  altitude.  As  men¬ 
tioned  previously,  neglecting  the  lateral  centrifugal  term  in  the  equa¬ 
tions  of  motion  may  result  in  a  higher  turn- angle  prediction.  This 

(13) 

error  can  be  30  percent,  as  pointed  out  by  Wang  for  the  severe  con¬ 
ditions  he  employed. 

Equation  (24)  can  be  rewritten  as 


gYJB (3  -  2) 


[vVF  -  73  -  2][VeV3~+yp  ‘  2] 
ivjs  +  Jr-  jar?!] 


(24a) 


But  </3(3  -  2)  3  for  large  values  of  3,  hence 


-  2)  » 


Lt  « 

u  —  cos  * 
o  D 


(24b) 


The  flight  times  for  various  flight  conditions  as  given  by  Eq.  (24) 
are  plotted  in  Fig.  6.  Substituting  Eq.  (24b)  into  Eq.  (24a)  shows 
that  t  is  proportional  to  L/D  cos  $.  In  other  words,  the  time  de¬ 
creases  to  zero  as  the  bank  angle  approaches  90  deg  or  L/D  approaches 
zero.  It  is  evident  that  high-entry-speed  and  high -L/D  vehicles 
require  a  longer  time  to  slow  down  to  the  same  speed  ratio.  For  a 
bank  angle  equal  to  90  deg,  the  equations  are  not  applicable,  since 
the  flight  path  would  resemble  a  ballistic  trajectory  rather  than  a 
lifting  trajectory. 

The  flight-path  length,  S,  for  two  different  flight-path  angles, 
Y  =  0.01  and  0.02,  is  shown  in  Figs.  7a  and  7b,  which  indicate  that 
more  than  half  of  the  distance  to  be  flown  will  be  in  the  region 
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where  the  flight  speed  is  greater  than  90  percent  of  the  entry  speed. 

In  general,  the  flight-path  length  increases  with  increasing  L/D  at 
a  given  speed  ratio. 

As  expected,  if  mean  value  y  is  used,  the  flight-path  length  is 
a  linear  function  of  altitude,  as  given  by  Eq.  (26)  and  shown  in  Fig. 

8.  It  can  be  seen  that  the  flight-path  length  would  be  doubled  when¬ 
ever  the  flight-path  angle  is  reduced  by  approximately  half. 

The  longitudinal  and  lateral  ranges  for  constant-bank-angle  flight 

are  given  by  Eqs.  (28)  or  (29)  and  (31)  or  (32),  respectively.  The 

numerical  results  for  L/D  =  1,  2,  3,  and  4,  and  V  =0.98  and  0.90, 

e 

are  presented  in  Figs.  9a  through  9e.  It  is  clear  that  large  ranges 
can  be  achieved  by  increasing  L/D.  The  maximum  lateral  range  can  be 
obtained  with  a  bank  angle  of  about  43  deg  at  L/D  =  1  and  V£  *  0.98. 
However,  this  bank  angle  for  maximum  lateral  range  shifts  to  about  33 
deg  at  L/D  =  3  at  low  speed.  The  vehicle  with  high  L/D  (>  3)  will  pro¬ 
ceed  along  its  spiralling  course  far  enough  to  achieve  large  heading 
changes  and  will  be  capable  of  reversing  its  direction  of  flight  as 
it  reaches  very  low  speed.  At  a  given  speed  ratio,  both  longitudinal 
and  lateral  ranges  are  very  sensitive  to  the  initial  entry  speed.  A 
comparison  of  Figs.  9c  and  9e  shows  that  for  y  =  0.01  and  L/D  =  3,  a 
reduction  of  entry  speed  from  0.98  uq  to  0.90  uq  may  reduce  the  longi¬ 
tudinal  range  by  more  than  half  and  may  cut  the  lateral,  range  by  more 
than  one-third  by  the  time  the  speed  ratio  of  0.1  is  reached. 

A  comparison  of  equilibrium-glide  results  obtained  in  Ref.  12 
(based  on  the  assumptions  of  a  very  small  flight-path  angle)  and  those 
of  the  present  analysis  (based  on  the  small  and  slowly  changing  flight- 
path  angle  of  Y  =  0.01  rad)  is  given  in  Table  1  and  Fig.  10.  These 
show  that  the  assumption  of  a  vary  small  flight-path  angle  results  in 
a  smaller  range  prediction  particularly  for  ;he  lateral  range  of  high- 
L/D  vehicles,  and  in  prediction  of  a  larger  bank-angle  requirement  for 
maximum  lateral  range- 

One  of  the  results  of  this  study  is  a  closed-form  solution  for 
predicting  the  performance  of  equilibrium-glide  vehicles  more  accu¬ 
rately  than  was  possible  in  the  past  without  machine  programming  of 
the  equations  of  motion  directly.  It  can  be  concluded  that  for 
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ToKl o  1 

X  V>  VXS.  A. 


EFFECT  OF  L/D  AND  FLIGHT-PATH  ANGLE  ON 
LONGITUDINAL  AND  LATERAL  RANGES3 
(V  =  0.98,  V  =  0.2V  ) 

2  O 


y 

yY*o.oi 

x  at 

y 

'max 
(n  mi) 

XY^-0.01 

$  at 

y 

'max 

(deg) 

SteO.Ol 

L/D 

(rad) 

max 
(n  mi) 

yY«0 

XY*=0 

$ 

Y«0 

i. 

0.01 

«;0 

580 

543 

1.068 

4,197 

3,750 

1.119 

42 

45 

0.933 

2 

0.01 

«i0 

2,236 

1,910 

1.171 

9,080 

7,650 

1.187 

38 

41 

0.927 

3 

0.01 

«s0 

4,780 

3,700 

1.292 

14,450 

11,200 

1.290 

33 

37 

0.892 

4 

0.01 

0 

8,940 

7,000 

1.277 

28,550 

15,700 

1.818 

27 

31 

0.871 

aBased  on  assumptions  of  a  very  small  flight-path  angle  (Y  «  0)  and 
small  and  slowly  changing  flight-path  angle  (Y  0.01  rad).  Data  for 
Y  R=!  0  were  obtained  from  Ref.  12. 


orbital-speed  reentry,  a  vehicle  with  L/D  =  3,  gliding  at  a  flight- 
path  angle  of  0.01  rad  and  a  constant  batik  angle  of  about  33  deg,  can 
provide  a  reasonable  lateral  range  greater  than  the  earth’s  radius. 

In  addition,  more  than  half  the  longitudinal  range  will  be  achieved 
before  the  flight  speed  reaches  90  percent  of  the  reentry  velocity. 
However,  half  the  lateral  range  will  be  achieved  when  the  vehicle 
velocity  is  approximately  60  percent  of  the  reentry  velocity  for  L/D 
=  1,  and  70  percent  for  L/D  =  3,  at  optimum  bank  angle. 
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